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Abstract 

The effects of combined external electric and magnetic fields on elastic collisions in ultracold 
Li-Rb mixtures is studied using recently obtained, experimentally verified potentials. Our analysis 
provides both quantitative predictions for and a detailed physical interpretation of the phenomena 
arising from electric-field-induced interactions. It is shown that the electric field shifts the positions 
of intrinsic magnetic Feshbach resonances, generates copies of resonances previously restricted to 
a particular partial-wave collision to other partial wave channels, and splits Feshbach resonances 
into multiple resonances for states of non-zero angular momenta. It was recently observed that 
the magnetic dipole-dipole interaction can also lift the degeneracy of a ^)-wave state splitting the 
associated p-wave Feshbach resonance into two distinct resonances at different magnetic fields. Our 
work shows that the splitting of the resonances produced by an applied electric field is more than 
an order of magnitude larger. This new phenomenon offers a complementary way to produce and 
tune an anisotropic interaction and to study its effect on the many-body physics of heteronuclear 
atomic gases. 
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I. INTRODUCTION 



The discovery of magnetic field tunable Feshbach resonances (FRs) has led to many 
groundbreaking experiments in the field of ultracold atomic and molecular physics M, 
3||. Magnetic FRsprovide a powerful tool to control microscopic interactions in ultracold 
quantum gases 0,01, offer an extremely_sensitive probe of interatomic interaction potentials 



8|, and can be used to create ultracold molecules 



for collisions at ultracold temperatures 
by coherently linking ultracold atoms [9|, llOj. FRs arise due to coupling between a quasi- 
bound molecular state in a closed collision channel and the scattering wave function of the 
colliding atoms in an open channel. Because the quasi-bound states and the free atomic 
pair have, in general, different magnetic moments, both their absolute energy and relative 
energy difference can be tuned using an external magnetic field. When the energy of the 
quasi-bound state is degenerate with the energy of the free atomic pair, a resonant scattering 
process occurs, the s-wave scattering length diverges and both elastic and inelastic collisions 
are dramatically enhanced. Recent theoretical work has also demonstrated the possibility 
of inducing FRs in heteronuclear mixtures of atomic gases by applying a static electric field 
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12l |. The mechanism is based on the interaction of the instantaneous dipole moment 



of the heteronuclear collision complex with the external electric field. This interaction is 
distinct from and, for bi-alkali mixtures with a large electric dipole moment, much larger 
than that responsible for electric field control of ultra-cold collisions based on the electric 
polarization of colliding atoms and the resulting dipole-dipole interaction 

Q,y,Q- The 

interaction considered here, for heteronuclear collisions, couples collision states of different 
orbital angular momenta, and the coupling becomes very significant near a FR. This coupling 
gives rise to new s-wave scattering resonances induced by the presence of FRs in higher 
partial wave states and can shift the positions of the quasi-bound states resulting in a 
shift in the positions of the intrinsic magnetic FRs (i.e. those present in the absence of 
an external electric field). Moreover, the electric field can induce a strong anisotropy of 
ultracold scattering by exerting a torque on the collision complex of ultracold atoms. 

The use of combined electric and magnetic fields to control interatomic and intermolecu- 
lar interactions has several distinct advantages over using magnetic fields alone. It has been 
demonstrated that the combination of electric and magnetic fields may be used to control 
both the position and width of FRs independently even for homonuclear collisions, leading 
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to complete control over the character of ultracold collisions [15||. In addition, for the rele- 
vant field strengths considered here (< 100 kV/cm and < 1 kG), the electric fields can be 
varied much faster than the magnetic fields, in large part, because the corresponding field 
energy density is 10 times smaller. Electric fields induce anisotropic interactions and angle- 
dependent scattering at ultra-cold temperatures, which may affect the dynamics of quantum 
degenerate gases in unanticipated ways. It is also interesting to note that the electric-field 
control of heteronuclear collisions can be achieved at fields low enough that they do not 
perturb the separated atoms or non-polar molecules since they interact significantly with 
electric field only when in a collision complex. For these reasons, electric field control of 
interatomic interactions may be preferable to magnetically or optically tunable scattering 
resonances in certain applications. 

Silber et al. have recently created a quantum degenerate Bose-Fermi mixture of ^Li and 
^^Rb atoms in a magnetic trap {l^ and detected inter-species FRs [l^l- These resonances 
may eventually be useful to improve the efficiency of sympathetic cooling in this mixture 
and for the study of boson-mediated BCS pairing jl^. These FRs also provide a method 
to create loosely bound LiRb dimers. LiRb molecules have a relatively large electric dipole 
moment (up to 4.2 Debye) lol, which makes the Li-Rb system a good candidate for research 
on ultracold dipolar gases and the experimental study of electric-field-induced FRs 0.Q. 
Recently, we have generated accurate singlet and triplet interaction potentials for ultracold 
collisions in ^Li-^''Rb gaseous mixtures by fitting the experimentally measured FRs In 
this work, we use these potentials to investigate in detail the effects of combined external 
electric and magnetic fields on elastic collisions in ultracold Li-Rb mixtures. To guide future 
experimental studies, we predict the positions and widths of electric-field-induced FRs for 
several spin states and explore the effect of the orientation of the electric field with respect 
to the magnetic field on ultracold elastic collisions. The work presented here represents the 
first quantitative analysis of electric-field-induced resonances based on precise inter-atomic 
potentials. In addition, our analysis provides insights into the detailed physical mechanism 
of electric-field-induced interactions in ultracold binary mixtures of alkali metal atoms. We 
report for the first time the observation that the coupling induced by electric fields splits FRs 
into multiple resonances for states of non-zero angular momenta. Recently it was observed 
that the magnetic dipole-dipole interaction can lift the degeneracy of a p-wave state 



This splits a p-wave FR into two distinct resonances at different magnetic fields 
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The splitting of the resonances studied here is produced only in heteronuclear collisions by 
a coupling between different quasi-bound states, is continuously tunable using an applied 
electric field, and is more than an order of magnitude larger than the splitting induced by 
magnetic dipole-dipole interactions. This new phenomenon offers a complementary way to 
produce and tune an anisotropic interaction and to study its effect on the many-body physics 
of heteronuclear atomic gases. 



II. THEORY 

The Hamiltonian for the ^Li-^^Rb collision system (or any bi-alkali) in the presence of 
superimposed electric and magnetic fields is 

H = 4ei + Ve{R) + Vb + (1) 

where ^rei accounts for the relative motion of the atoms, Ve{R) describes the interaction 
between LiRb and an external electric field, Vb models the interaction of the collision com- 
plex with an external magnetic field, and Vhf represents the hyperfine interactions. Hre\ can 
be written as 

- -2ilRdR^^ + ^lHW + ^^^^ 

where /z is the reduced mass of the colliding atoms, R is the interatomic distance, I is the 
operator describing the rotation of the collision complex and the angles 9 and (j) specify the 
orientation of the interatomic axis in the space-fixed coordinate frame. Here, we neglect the 
magnetic dipole-dipole interaction since it has a negligible effect on the observables described 
in this paper. The atomic and molecular quantum numbers used in this article are defined 
in Table [D 

We expand the total wave function of the diatomic system in a fully uncoupled, space- 
fixed basis set: 

a I mi 

where FaimiiR) \lmi) and |a) = \IuMj^.)\SuMs^,)\lRhMj^J\SKhMs^J are the radial basis 
states and the atomic spin states, respectively. The substitution of this expansion in the 
Schrodinger equation with the Hamiltonian (1) results in a system of coupled differential 
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equations 

(P l{l + 1) 



2/xe 



Falmi (-R) 



2^iJ2J2J2{almi V{R) + Ve{R) + Vj, + V^, a'l'm[) F^,vm[{R) (4) 



which we solve at fixed values of total energy e. 

The electronic interaction potential V{R) can be represented as 

V{R) = J2Y1 \SMs)Vs{R){SMs\ (5) 

S Ms 

where Vs{R) denotes the adiabatic interaction potential of the molecule in the spin state S. 
To evaluate the matrix elements of the interaction potential V{R), we write the atomic spin 
states \IuMI^.)\SuMs^■)\IKbMI^^^)\S1ihMsJ^^) in terms of the total electronic spin 

\IuMjJ\SuMsJ\IkM^,)\Ski.Ms,,) = 

5^5^(-l)^^^(25+l)i/M "^"^ "^^^ ^ ]\IuMjJ\I^M,,)\SMs) (6) 

\ Ms,, Ms,, -Ms J 

and note that 

{SMs\ViR)\S'M's) = Vs{R) Sss'hisM',- (7) 

The term enclosed in parentheses in Eq. [6] denotes a 3j-symbol. The operator V^(-R) is 
diagonal in the nuclear spin states and / and nii quantum numbers. 
The operator Ve{R) can be written in the form 

Ve{R) = -E ■ d= -E{eE ■ed)J2Yl \SMs)ds{R){SMs\ (8) 

S Ms 

where ce and id are the unit vectors pointing along the electric field and dipole moment 
of the LiRb dimer respectively, ds denotes the spin dependent dipole moment functions of 
LiRb, and E is the electric field magnitude. Clearly, the electric field coupling depends on 
the orientation of the field with respect to the dipole moment. Specifically, ■ = cos(x) 
where x is the angle between E and d. If the electric field and dipole moment vectors are 
oriented at angles 7 and 9 with respect to the quantization axis (taken to be along the 
i-axis), then this term can be written in terms of the first-degree Legendre polynomial as 

eE-ed = cos(x) = Pi(cos(x)) (9) 



where the l^x are spherical harmonics and we have used the spherical harmonic addition 
theorem. The azimuthal angles, (p^ and 06i, are measured from the positive x-axis to the 
orthogonal projection of the E and d vectors in the x-y plane. Figure [T] illustrates this 
coordinate system. 

The dipole moment functions are modeled by 

ds{R) = D exp [-a(R - Re)^] (10) 

with the parameters Re = 7.2 ob, a = 0.06 a^^ and D = 4.57 Debye for the singlet state, 
and Re = 5.0 ob, a = 0.045 a^^ and D = 1.02 Debye for the triplet state, where the Bohr 
radius is ob = 0.0529177 nm. These analytical expressions approximate the true functions 
are were fit to the numerical data for the dipole moment functions computed by Aymar and 
Dulieu 0. 

The collision dynamics in Li-Rb system depends on the relative angle (7) between the 
electric and magnetic fields. Without loss of generality, we can assume that E lies in the 
x-z plane, i.e., (p^ = 0. The matrix elements of Ve(-R) are therefore evaluated using the 
expressions 



{lmi\eE ■ ed\l'm[) ■ 

v2 



1 , , I I I l'\ I I 1 I' \ 

-=sin7(-irV(2/ + l)(2/' + l) (11) 
/2 \ / \ mz -1 -m[ I 



+ cos7(-l)'"V(2^+l)(2/' + l) 



i=sin7(-l)™V(2/ + l)(2/' + l) 




and 



{SMs\ I J2 Yl \S"M'^)ds"{S"M'^\ I \S'M's) = ds 6ss'SmsM',. (12) 



S" M'^ 



It is important to note that the electric field coupling operator, Ve{R), has a spin structure 
identical to the electronic interaction potential, V{R). Namely, they are both diagonal in 
the total electronic spin and its z projection. However, as is clear from Eq. [H], the geometric 
factor, eE-ed, introduces this additional electric-field-induced coupling only between states of 
different orbital angular momenta. The first 3j-symbol in Eq. [TT]is non-zero for / + /' = ±1. 
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The second 3j-symbol in Eq. [TT] selects states with orbital angular momenta projections 
differing by Ami = m/ — mj = ±1, 0. If the electric field is directed along the 2;-axis (the 
quantization axis) then 7 = and Eq. [TT] reduces to 

, / / 1 /' W / 1 /' \ 

{lmi\ cose\l'm'i) = 5„,^;(-l)™V(2^ + 1)(2/' + 1) . (13) 

y -mi mi J \^000y 

In this case only those states with the same orbital angular momenta projections are coupled. 

In the absence of an electric field, the Hamiltonian (neglecting magnetic dipole-dipole 
interactions) preserves the projection of the total angular momentum along the quantiza- 
tion (magnetic field) axis, i.e. the sum Ms^- + Ms^^ + Mj^^ + Mj^^ + mi is conserved. In 
addition, the orbital angular momentum is also conserved. Therefore, the intrinsic {E = 0) 
FRs which are present for a particular atomic spin state arise from its coupling to the set 
of bound states with the same total angular momentum projection and orbital angular mo- 
mentum. The energy of the bound states relative to the threshold state and their relative 
magnetic moments determine the magnetic fields at which resonances will occur and the 
strength of the coupling to a particular bound state determines the width of the resulting 
resonance. The electric field induces additional couplings between the threshold state and 
bound states of different angular momenta but still within the same set of states sharing the 
same total spin angular momentum projection mp = M^^; +Ms^^ +^iLi + ^^Rb- The result 
is that FRs previously restricted to a particular partial- wave collision will appear on adja- 
cent partial wave scattering states. The widths of these electric-field-induced FRs depend 
on the strengths of the couplings and therefore on the magnitude of the electric field, E. In 
addition, the electric field coupling among the bound states (with the same total mp value) 
will shift their absolute energies and therefore the positions of the intrinsic FRs. Moreover, 
this coupling and the resulting shifts are, in general, dependent on mi and, as a consequence, 
resonances associated with bound states with / > will split into / + 1 distinct resonances. 

The interaction of the atoms with an external magnetic field B is described by 

Vs = 2/io5 + 5zh.) - B (iJ/z,. + ^4,,) (14) 

where B is the magnetic field strength (directed along the z-axis), //q is the Bohr magneton 
and /iLi(Rb) denote the nuclear magnetic moments of Li (Rb) . Sz^.^^^^^ and Izi^^j^^ give the 
z component of the operators describing the electronic and nuclear spins of Li (Rb), 5'Li(Rb) 
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and /Li(Rb), respectively. The hyperfine interaction Vhf can be represented as 

^hf = luJu ■ Su + TRb^Rb ■ 'S'Rb (15) 

where and 7Rb are the atomic hyperfine interaction constants: = 152.14 MHz and 
7Rb = 3417.34 MHz (where we work in units with h = 1). 

The operator representing the magnetic field interaction is diagonal in the representation 
|/LiM/j^.)|5'LiM5'Li)|/Rb^7Hb)l'S'Rb^5Rb)i while the matrix elements of the hyperfine interac- 
tion operators are not but can be readily evaluated using the relations: 

-^Li ■ 'S'Li = /zli'S'zl; + -{iu+Su- + Iu-Su+) (16) 

and 

-^Rb ■ S^h = -^^Rb'S'zRb + 2 (-^Rb+'^R-b- + -^Rb-'S'Rb+)- (17) 

I± and S-i- are the raising and lowering operators. Because of these terms, the matrix of the 
Hamiltonian in the basis \almi) does not become diagonal as i? ^ oo. Therefore, we find 
the basis which diagonalizes the matrix of Ve + Vb + Vhf, and transform the solutions of 
Eq. [4] into this basis. At this point, the boundary conditions are applied and we construct 
the scattering S'-matrix. This procedure has been described in Ref. [22||. The scattering 



matrix thus obtained yields the probabilities of elastic and inelastic scattering of Li and Rb 
in the presence of electric and magnetic fields. 

III. RESULTS 

The Li-Rb mixture is an important system for the study of both ultracold atomic and 
molecular gases. A quantum degenerate Bose-Fermi mixture of ^Li and ^^Rb atoms has 

n . 

been recently created |16l | and may be important for research of ultracold fermionic and 
bosonic mixtures. By tuning an external magnetic field through a FR, researchers can 
create an ensemble of LiRb molecules, a polar dimer, from this gaseous mixture. Since LiRb 
molecules have a large electric dipole moment, the Li-Rb system is also a good candidate for 
the research of ultracold dipolar gases and the experimental study of electric-field-induced 



FRs [12||. Motivated by these features and the availability of accurate potentials [7||, we 
analyze the collision properties of Li-Rb mixtures in the presence of both magnetic and 
electric fields. 
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As discussed above, in the absence of electric fields, diff'erent partial wave states \lmi) 
of the Li-Rb collision complex are uncoupled and s-wave scattering entirely determines 
the collision dynamics in ultracold Li-Rb gases. The presence of an external electric field, 
however, induces couplings between states of different orbital angular momenta with Al = 
±1. As a result, a resonant enhancement of the s-wave cross section appears at magnetic 
fields near intrinsic p-wave resonances. Figure [2] shows the magnetic field dependence of s- 
and p-wave elastic cross sections for Li and Rb atoms in the spin state |^, ^)6u ® |1) l)87Rb 
(where \f,mf) is the usual notation for the atomic hyperfine states) computed at zero electric 
field and at E = 100 kV/cm. Here the electric field is directed along the quantization axis 
(7 = 0). In the presence of the 100 kV/cm electric field, an s-wave resonant peak appears 
(indicated at A) at the magnetic field of 877.50 G arising from an intrinsic p-wave resonance. 
We refer to this resonance as an electric-field-induced Feshbach resonance. Figure [3] shows 
this feature in more detail. We present in Table HTl the positions and widths of electric-field- 
induced FRs for several atomic spin states of the ®Li-^^Rb system at magnetic fields below 
2 kG. For each resonance, we extract the position (Bq) and width (AB) from the magnetic 
field dependence of the scattering length. 

a(S) = abg (1 - -T^^] (18) 



where a^g represents the background scattering length. In this calculation, we also observe 
new p-wave resonances induced by the coupling to a d-wave state, and we find that these 
p-wave resonances give rise to new s-wave electric-field-induced FRs, denoted by (d) in 
Table M 

The width of the electric-field-induced FRs is determined by the strength of the coupling, 
which is, in turn determined by the magnitude of the electric field. In Fig. HJ we plot the 
width of the s-wave resonance (shown in Fig. [3]) induced by the intrinsic p-wave resonance 
near 882 G as a function of the electric field magnitude. We find that the width can be well 
represented by a quadratic function of at least for the electric fields below 200 kV/cm, 
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which suggests that this induced resonance arises from an indirect coupling 

The electric field not only induces new resonances but also shifts the position of intrinsic 
magnetic FRs. Figure [2] shows that the interaction of Li-Rb dipole moment with the electric 
field shifts the position of both the s- and p-wave resonances. At B and at C an intrinsic s- 
wave resonance is shifted to higher magnetic fields (corresponding to a shift of the associated 
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bound state to lower energy) due to the electric field. At D an intrinsic p-wave resonance 
is shifted to lower magnetic fields (corresponding to a shift of the associated bound state to 
higher energy). For the most part, the shift of the FR positions arises from the coupling 
between bound states whereas the coupling of a given bound state to the scattering state 
results in a broadening of the associated electric-field-induced resonance. We note that the 
shift of magnetic FRs at higher magnetic fields (e.g. C) is more significant than the shift at 
lower magnetic fields (e.g. B). This generic behavior results from the fact that resonances 
associated with higher magnetic fields are typically more deeply bound than those associated 
with lower magnetic fields. As a result, the wave function of the bound state giving rise to 
FRs at higher fields samples smaller interatomic distances where the dipole moment function 
is much larger. 

Another example of the shift induced by the electric field couplings is shown in Fig. [5] for 
atoms in the atomic spin state ||, — |)6Li ® |1, — l)87p^b. An electric field of 30 kV/cm is large 
enough to shift the position of this s-wave resonance by almost 2 G - much larger than its 
width - while a field of 100 kV/cm produces a shift of almost 9 G. It is important to note that 
this s-wave resonance shifts to lower magnetic fields as the electric field increases, and this 
is opposite to the shift of the s-wave resonances shown in Fig. [2l The shift of a resonance 
results from level repulsion between the closed channel bound states and therefore depends 
on the proximity, position, and coupling strengths of the nearby bound states. Therefore, 
the direction of the resonance shift and its dependence on the electric field magnitude do not 
exhibit a generic behavior but depend on the particular environment of a given resonance. 

These shifts provide a way to dramatically and rapidly modify the s-wave scattering 
length by tuning into and out of an intrinsic magnetic field resonance. Fig. [6] presents two 
such tuned resonances arising from the variation of the electric field for atoms in the atomic 
spin state ||, |)6Li ® |1, l)87j^b. This figure shows the cross section for s-wave collisions as 
a function of the electric field strength with the magnetic field fixed at 1066 G (solid line) 
and 878 G (dotted line). The solid curve shows a large resonance feature due the intrinsic 
magnetic FR at 1067 G which shifts to higher magnetic field (lower absolute energy) as 
the electric field increases (see Fig. [2|). The small resonance feature which appears at the 
electric field strength of approximately 16 kV/cm in the solid curve arises from an electric- 
field-induced resonance arising from the intrinsic p-wave resonance just above 1066 G which 
shifts to lower magnetic field (higher energy) as the electric field increases. In the same 
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plot, the dotted curve shows a resonance feature due to the shift of an electric-field-induced 
resonance arising from the intrinsic p-wave resonance at 882 G. Fig. [3] shows the p-wave state 
responsible for this resonance shifts to lower magnetic fields (higher energy) as the electric 
field increases. 

In the presence of an electric field, the couplings between different partial wave states 
can push the bound states in, for example, the s-wave and p-wave interaction potentials 
apart. This level repulsion gives rise to the electric-field-induced shift of the intrinsic s- 
and p-wave magnetic FRs. Since the couplings depend on the orbital angular momentum 
projection, rrii, we also expect the electric field induced coupling to also split FRs for states 
of non-zero angular momenta. This mechanism is illustrated in Fig. [7| where three adjacent 
bound state levels are shown as well as the coupling induced by an applied electric field 
with 7 = 0. Without external electric fields, the bound states in the p-wave interaction 
potential are degenerate, whereas the electric field lifts this degeneracy. In the case where 
the electric field points along the quantization axis (7 = 0), the rrii = bound state in the 
p-wave potential is coupled to bound states in both the s- and rf-wave potentials whereas the 
\mi\ = 1 bound states are only coupled to bound states in the (i-wave potential. This occurs 
because the system is cylindrically symmetric and the couplings between internal spin states 
and partial wave states are negligible. As a result, the coupled states repel and the rrii = 
state is shifted differently than the |m/| = 1 states splitting the p-wave resonance into a 
doublet. For the purposes of simplifying the discussion, we have neglected the coupling to 
yet higher order partial wave states and we have neglected the possible presence of other 
closed channel states in the near vicinity. This mechanism generally applies to all nonzero 
partial waves. For a state with an orbital angular momentum /, the number of peaks is 
/ + 1 corresponding to the number of distinct values for |m/|. Fig. [8] and Fig. [9] show the 
splitting of a p-wave and a d-wave FR, respectively. In the presence of a 100 kV/cm electric 
field, the p-wave resonance splits into two peaks (corresponding to the |m;| = 1 and m; = 
components) with a separation of 4 G (dash-dot line in Fig. [8]). The shift of the \mi\ = 1 
peak in Fig. [8] to higher magnetic fields (lower energy) is consistent with coupling between 
the p-wave bound states and a d-wave state which resides at a higher energy (illustrated 
in Fig. ini). The splitting of a d-wave bound state gives rise to three separated resonances 
and is shown in Fig. [9l An interval of 1 G opens up between = and \mi\ = 1 and 
the an interval of 2 G appears between \mi\ = 1 and \mi\ = 2. Since it is only very weakly 
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coupled to higher partial- wave states, the \mi\ = 2 component remains in essentially the 
same location as the resonance at zero electric field. 

The splitting of FRs in states with non-zero orbital angular momenta has been previously 
discovered in experiments by Jin and coworkers [2^ and studied theoretically by Ticknor et 
al. I2II. They found that p-wave FRs for collisions of homonuclear gases of split into 
a doublet due to the magnetic dipole-dipole interaction. In the work presented here, we 
neglect the magnetic dipole-dipole interaction since it produces a negligible effect compared 
to the electric field coupling and the splitting we predict for FRs is entirely due to the effect 
of the electric field. As discussed in 2^, the ability to introduce and tune an anisotropic 
interaction using Mgh-pa.tia.-wave resonances may have ta. paching consequences to. the 
study of novel forms of superfiuidity using cold atomic gases [24] . 

The splitting of the nonzero-partial- wave resonances arising from magnetic dipole-dipole 
interactions is very small and will disappear as the resonance becomes broad with increasing 
temperature. In contrast, the splitting observed here, occurring for heteronuclear atomic 
mixtures, is more than an order of magnitude larger. This new phenomenon offers a com- 
plementary way to produce and tune an anisotropic interaction and to study its effect on 
the many-body physics of heteronuclear atomic gases. 

So far, we have discussed the modifications of FRs induced by the application of an 
electric field parallel to the magnetic field (7 = 0). In addition, we study the effect of non- 
parallel fields (7 7^ 0). In Fig. [TOl we show the variation of the total elastic cross section for 
p-wave collisions given fixed electric (100 kV/cm) and magnetic fields as a function of the 
angle between them 7. In the upper panel, the magnetic field was 877 G which is near the 
p-wave resonance for the m; = component (see Fig. [8|). Whereas, in the lower panel, the 
magnetic field was 881.9 G and falls in between the m/ = and \mi\ = 1 resonances in the 
p-wave doublet. In the latter case, the variation of the cross sections as a function of 7 is 
only a factor of 10, while at a magnetic field near one of the resonances, the cross section 
varies by almost 4 orders of magnitude as 7 changes by less than 30°. 

Fig.dUpresents the magnetic field dependence of the total elastic cross section for different 
components of p-wave scattering at 7 = 45° near the intrinsic p-wave resonance at 882 G. In 
this case, because the electric field couples states of differing rrii values, the doublet structure 
of the p-wave resonance appears on each of the three nii components of the open channel. 
This is in contrast to the case with 7 = shown in Fig. [8] where the coupling is only between 
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states with the same mi value and each component exhibits a single resonance. It should 
be clarified here that the electric-field-induced s-wave resonance arising from this p-wave 
resonance exhibits only the single resonance corresponding to the m; = component of 
the p-wave bound state. This is because (neglecting the magnetic dipole-dipole interaction) 
the orbital angular momentum projection along the electric-field axis is conserved by the 
Hamiltonian, and m; = for s-wave collisions in all coordinate frames. On the other hand, 
a state with orbital angular momentum / and projection mi defined with respect to the 
magnetic field axis will be a linear combination of states with all possible values of mi when 
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represented with respect to the electric field axis 

Fig. [I2] presents the magnetic field dependence of the average p-wave elastic scat- 
tering cross section (averaged over all three components) for atoms in the spin state 
||)|)i3Li ® |l)l)8''Rb at -E = 100 kV/cm and with three orientations of the electric field, 
7 = 0°, 45°, and 90°. The main point of this plot is to illustrate that the position of the 
resonances remains unchanged for different values of 7. This is particularly important for 
the experimental search for these effects since it means that any variation of the orientation 
of the electric and magnetic fields does not adversely affect the visibility of these multiplet 
features. Consequently, any inhomogeneities in the direction of the electric field over the 
confinement size of the atomic ensemble would also not affect their visibility. Of course, 
since the positions of the resonances do depend on the electric field strength, any inhomo- 
geneities in the magnitude of the electric field would result in inhomogenous broadening of 
the observed resonances. 



IV. CONCLUSIONS 



We have presented quantitative predictions of the eff'ects of combined external electric 
and magnetic fields on elastic collisions in ultracold Li-Rb mixtures. This work is the first 
analysis of electric-field-induced interactions based on precise, experimentally verified inter- 
atomic potentials. In addition, we have provided important insights into the detailed physical 
mechanism of electric-field-induced interactions in ultracold binary mixtures of alkali metal 
atoms. We have shown that the electric field both shifts the position of intrinsic FRs and 
generates copies of resonances previously restricted to a particular partial-wave collision to 
other partial wave channels. To facilitate the experimental search for these phenomena, 
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we have provided predictions for the positions and widths of electric-field-induced FRs for 
several spin states and we have analyzed the effects of a non-parallel orientation of the electric 
field with respect to the magnetic field on ultracold elastic collisions. We also have reported 
for the first time the observation that the coupling induced by electric fields splits FRs into 
multiple resonances for states of non-zero angular momenta. It was recently observed that 
the magnetic dipole-dipole interaction can also lift the degeneracy of a p-wave state splitting 



the associated p-wave FR into two distinct resonances at different magnetic fields [20|, l2l| . 
The primary differences with that work are that the splitting studied here is produced only 
in heteronuclear collisions, is continuously tunable using an applied electric field, and is 
more than an order of magnitude larger than the splitting induced by magnetic dipole- 
dipole interactions. We believe the additional degrees of control offered by electric-field 
interactions will play an important role in future experiments on the many-body physics of 
heteronuclear atomic gases. 
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TABLE I: Definition of quantum numbers used in this paper. 



/ 


orbital angular momentum of the diatomic system 


mi 


projection of I on the space-fixed quantization axis 


S 


total electronic spin angular momentum of the diatomic system 


Ms 


projection of S on the space-fixed quantization axis 




nuclear spin angular momentum of Li 




projection of lu on the space-fixed quantization axis 


Im, 


nuclear spin angular momentum of Rb 




projection of on the space-fixed quantization axis 


'S'Li 


electronic spin angular momentum of Li 


Ms,, 


projection of S'Li on the space-fixed quantization axis 




electronic spin angular momentum of Rb 




projection of on the space-fixed quantization axis 
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TABLE II: The positions (Bq) and widths (Ai?) of s-wave resonances induced by an external electric 
field of 100 kV/cm for ^Li-^^Rb at magnetic fields below 2 kG. {d) denotes an s-wave electric-field- 
induced Feshbach resonance arising from a second order coupling through the p-wave channel to a 
d-wave closed channel state. As a consequence, these resonances are exceedingly narrow. 





Atomic States 


Bo 


AB 


\f,'mf)6u 


1/, m/)87Rb 


(G) 


(G) 


l2' 2/ 




536.65 (d) 
877.5 

654.52 


0.01 
2.3 
< 0.01 


l2' 2/ 


|i,o) 


555.88 {d) 
885.8 


< 0.01 
2.6 


l2' 2/ 




578.58 (d) 
707.70 
770.50 


0.01 

< 0.01 

< 0.01 


li _1\ 
l2' 2/ 


|1>0) 


596.01 (d) 
926.8 


< 0.01 
2.6 


|3 3\ 
l2' 2/ 




1242.5 


12.7 
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X 



FIG. 1: (Color online) The coordinate system used for our calculations. E and d* represent the 
vector of the external electric field and the dipole moment vector, respectively; 7 specifies the angle 
between the electric field vector and the quantization axis; 9 is the angle between the dipole moment 
vector and the z-axis; and x is the angle between E and d. The azimuthal angles, cj)^ and 4)q, are 
measured from the positive x-axis to the orthogonal projection of the E and d vectors in the x-y 
plane. 
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900 1000 1100 1200 1300 

Magnetic field (G) 

FIG. 2: (Color online) Magnetic field dependence of the elastic cross section for collisions between 
Li and Rb in the atomic spin state ^)6Li ® |1, l)87j^b. These results were obtained for a collision 
energy of 10~''cm~^ and two different electric fields. The solid and dash-dotted curves show the s- 
and p-wave cross sections with E = 0, while the dotted and dashed curves show the s- and p-wave 
cross sections when E = 100 kV/cm. Here, only the cross section for the m; = state is shown for 
p-wave scattering. At A an s-wave resonance is induced by an intrinsic p-wave resonance. Figure 
[3] shows this feature in more detail. At B and at C an intrinsic s-wave resonance is shifted to 
higher magnetic fields (corresponding to a shift of the associated bound state to lower energy) due 
to the electric field coupling between bound states. The observation that the shift of higher field 
resonances (e.g. C) is typically larger than that of lower field resonances (e.g. B) is discussed in 
the text. At D an intrinsic p-wave resonance is shifted to lower magnetic fields (corresponding to a 
shift of the associated bound state to higher energy). At E an induced p-wave resonance appears 
(invisible on this scale) due to the intrinsic s-wave resonance at C. 
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Magnetic field (G) 



882 



FIG. 3: (Color online) Magnetic field dependence of s- and p-wave elastic cross sections for atoms 
in the atomic spin state ^)6Li |1, 1)8?^^ computed at different electric fields. This is the same 
feature at A in Fig. [2l The solid and dotted curves show the s-wave cross sections at E = and 
E = 100 kV/cm, respectively. The dot-dashed and dashed curves show the p-wave cross sections at 
E = and E = 100 kV/cm, respectively. This intrinsic p-wave resonance shifts to lower magnetic 
field (corresponding to the shift of the associated bound state to higher energy) as the electric field 
magnitude is increased. The s-wave induced resonance appears at the same location as the intrinsic 
p-wave resonance, and its width grows with the strength of the electric field (see Fig.[l|. Here only 
the cross section of the = component is shown for the p-wave state is shown (Fig. [8] shows the 
cross sections for all three components). The collision energy is 10~^cm~^. 
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250 



Electric field (kV/cm) 



FIG. 4: (Color online) The width (AB) of the s-wave electric-field-induced Feshbach resonance 
arising from the intrinsic p-wave resonance at 882 G as a function of the electric field magnitude. 
Here 7 = and the collision energy is 10~^cm~^. The width appears to scale quadratically with 
at least for the electric fields below 200 kV/cm, and suggests that this induced resonance arises 



from an indirect coupling 
of kV/cm. 



2311. The solid line is the fit A5 = 1.76 x 10""^ E"^ G, where E is in units 
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FIG. 5: (Color online) Magnetic field dependence of the s-wave elastic cross section for atoms in 
the atomic spin state |^,— ^)6Li (E) |1,— l)87j^b computed at different electric fields: E = kV/cm 
(solid curve), = 30 kV/cm (dotted curve), E = 70 kV/cm (dashed curve) and E = 100 kV/cm 
(dot-dashed curve). An intrinsic s-wave resonance (whose position is 1611 G in the absence of an 
electric field) is observed to shift to lower magnetic fields as the electric field strength is increased. 
Note: the shift direction is in the opposite sense to that of the intrinsic s-wave resonances in Fig. El 
These results were obtained with a collision energy of 10~''cm~^. 
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Electric field (kV/cm) 

FIG. 6: (Color online) Variation of the cross section for s-wave collisions as a function of the electric 
field strength with the magnetic field fixed at 1066 G (solid line) and 878 G (dotted line) for atoms 
in the spin state ^)6u^ |1, l)87j^b- The large resonance feature shown in the solid curve is due to 
the shift of the intrinsic magnetic Feshbach resonance just below 1066 G to higher magnetic fields, 
while the small resonance feature at 16 kV/cm arises from the shift of an intrinsic p-wave resonance 
just above 1066 G to lower magnetic fields as the electric field increases. The dotted curve shows a 
resonance feature associated with an electric-field-induced resonance (shown in Fig. [3]) which moves 
from 882 G at = down to a magnetic field below 877 G at -B = 120 kV/cm. The collision 
energy is 10" ''cm"^. 
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FIG. 7: (Color online) A schematic illustrating the mechanism of the splitting of p- and d-wave 
bound states resulting in the splitting of the corresponding Feshbach resonances. For simplicity, 
only three adjacent bound state levels are shown. The different partial wave potentials of each 
state are on this scale almost indistinguishable and are drawn here as a single potential. The inset 
shows the energy levels associated with these three states. The dotted lines indicate their energies 
in the absence of an electric field. The coupling induced by the electric field is represented as 
double-ended arrows and shown for the case when the electric field is aligned along the magnetic 
field, i.e. when 7 = states with the same m; value are coupled. The coupling results in level 
repulsion and the new position of the states is indicated by the solid lines. The degeneracy of the 
p- and d-wave bound states is broken and the associated Feshbach resonance splits into a multiplet 
with I + 1 distinct resonances as shown in Figs. [8] and [9l This simple picture predicts that the 
s-wave resonance should shift to higher magnetic fields (given the energy of the threshold moves 
down with increasing magnetic fields) and that the m; = partial wave component should produce 
a new resonance at a magnetic field below the^; = 1 component - consistent with the motion of 
the resonances in Fig. [2] and Fig.[8l Of course, each state is coupled to all other bound states within 
the same spin manifold and with an orbital angular momenta differing by A/ = ±1, resulting in 
splittings and shifts (e.g. Fig. [5]) which may not follow the predictions of this simple picture 




Magnetic field (G) 

FIG. 8: (Color online) Magnetic field dependence of p-wave elastic cross section (averaged over all 
three orbital angular momentum components) for atoms in the atomic spin state ^)6LiC5i|l, l)87Rb 
computed at zero electric field (solid curve) and at E = lOOkV/cm (dot-dashed curve). The thin 
dotted curves show the magnetic field dependence of the cross section for the |m/| = 1 and the m = 
components separately. The p-wave resonance splits into two distinct resonances, one occurring for 
the m/ = component and one for the |m/| = 1 components. When the electric and magnetic fields 
are not co-linear, this segregation of the resonance multiplet breaks down as seen in Fig. [T21 The 
collision energy is 10~^cm^^. 
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: E = OkV/cm 




Magnetic field (G) 

FIG. 9: (Color online) The upper panel shows the magnetic field dependence of the d-wave elastic 
cross section for atoms in the atomic spin state ^)6u ® |1) l)**''Rb computed at zero electric fields 
(dotted-dashed curve). The lower panel shows the magnetic field dependence of d-wave elastic cross 
section (solid curve). The contributions to the cross section from the |m;| = 2, Im^l = 1 and the 
m = components are shown (dotted curves) at = 100 kV/cm. The d-wave resonance splits into 
/ + 1 = 3 distinct resonances corresponding to the splitting of the d-wave bound state levels drawn 
schematically in the lower panel. The collision energy is lO^^cm^^. 
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FIG. 10: (Color online) Total elastic cross section for different components of p-wave scattering 
versus the angle, 7, between the applied electric and magnetic fields. The cross sections are shown 
for collisions in the m/ = state (dashed curve), the \mi\ = 1 states (dotted curve), and the average 
(solid curve) of the cross sections over all three components for the atomic state |)6Li(Xi|l, 1)87^^ 
and for E = 100 kV/cm . The upper panel shows these cross sections at an applied magnetic field 
of 877.0 G which is near the resonance for the m; = component while the lower panel is at a field 
of 881.9 G which is in between the resonances for the m; = and |m/| = 1 components (see Fig.[8|). 
We observe that the shape of this variation changes dramatically near a resonance. The collision 
energy is 10~^cm~^. 
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Magnetic field (G) 

FIG. 11: (Color online) Magnetic field dependence of the elastic cross section for different compo- 
nents of wave scattering with an electric field, E = 100 kV/cm, titled with respect to the magnetic 
field axis by 7 = 45°. The cross sections are shown for collisions in the mi = state (dashed curve), 
the |m;| = 1 states (dotted curve), and the average (solid curve) of the cross sections over all three 
components for the atomic state ^)6Li'X> |1, 1)87^^- The doublet structure of the p-wave resonance 
seen also in Fig. [8] now appears on each of the three angular momentum projection components. 
The collision energy is 10~^cm~^. 
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FIG. 12: (Color online) Magnetic field dependence of elastic cross sections for atoms in the atomic 
spin state |^,^)6Li (X" |l,l)87j^^3 computed at E = 100 kV/cm with the orientation of the electric 
field at 7 = 0° (solid curve), 45° (dotted curve), and 90° (dot-dashed curve). The collision energy 
is 10~''cm~^. 
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